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Abstract 

I review recent progress on bottom-quark fragmentation in e + e~ annihilation, top-quark decay and 
, H — > bb processes. I discuss the implementation of collinear and soft resummation, and the inclusion 

■ of non-perturbative information taken from LEP and SLD data. 

Heavy-quark physics is currently one of the major fields of investigation in experimental and theoretical 
particle physics. In order to investigate heavy-quark phenomenology, fixed-order calculations are reliable 
as long as one considers inclusive observables, such as total cross sections. For less inclusive quantities, 
one needs to resum large contributions, which correspond to collinear or soft parton radiation. 

In this talk I will investigate bottom-quark production in e + e — > bb processes, top-quark decay 
t — > bW, and Standard Model Higgs decay H — ► bb. In particular, I wish to study the energy distribution 
of b quarks and 6-flavoured hadrons in such processes. The results will be expressed in terms of the b 
\ normalized energy fraction Xb, which is given, e.g. in Higgs decay H — > bb{g), by: 2 
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The differential distribution for the production of a massive b quark at next-to-leading order (NLO) in 
the strong coupling constant a,s is given by the following general result: 
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In Eq. (2) Tq is the width of the Born process, /i is the renormalization scale, A(xb) is a function 
independent of nib, p > 1, P qq {xb) is the Altarelli-Parisi splitting function: 

P qq (x b ) = C F ( 1 -±4) ■ (3) 



Equation (2) is formally equal to the bottom energy spectrum in top decay or e + e~ processes: one will 
just have to replace m# with the top mass or the centre-of-mass energy y/s of the e + e~ collision. Function 
A(xb) is instead process-dependent. 

Equation (2) presents a large mass logarithm ~ asln(m^/mj), which needs to be resummed to 
all orders to improve the prediction. To achieve this goal, we can follow the approach of perturbative 
fragmentation functions [2], which expresses the energy spectrum of a heavy quark as the convolution 
of a coefficient function, describing the emission of a massless parton, and a perturbative fragmentation 
function D(nib, (J-f), associated with the transition of a massless parton into a massive quark: 
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In Eq. (4) , dTi / dz is the differential width for the production of a massless parton i in Higgs decay with 
an energy fraction z; Di(x, up, mb) is the perturbative fragmentation function for a parton i to fragment 

x Talk given at 18th International Workshop on High Energy Physics and Quantum Field Theory (QFTHEP 2004), 
St.Petersburg, Russia, 17-23 June 2004. 

2 In Eq. (1) powers of m b /m,fi are neglected. In top decay we shall have an extra term to account for the W mass [1], 
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into a massive b quark, \ip is the factorization scale. Neglecting g — > bb splitting, i = b on the right-hand 
side of Eq. (4). The coefficient function for H — > bb processes has been computed in [3] and reads, in the 
MS factorization scheme: 
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As pointed out in Ref. [3] , in the calculation of Eq. (5) , using the MS -renormalized Yukawa coupling yb(rn,b) 
for the Hbb vertex turned out to be essential. The MS NLO coefficient function for e + e~ annihilation 
was calculated in [2], for top decay in [1]. 

The perturbative fragmentation function follows the Dokshitzer-Gribov-Altarelli-Parisi (DGLAP) 
evolution equations [4, 5]. Its value at a given scale \xf can be obtained once an initial condition is given. 
In [2] the initial condition D™ 1 (xb, A*of, m-h) was calculated and its process-independence was established 
on more general grounds in [6]. It is given at NLO by: 
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The NNLO initial condition was calculated in Refs. [7, 8]. 

As discussed in [2], solving the DGLAP equations for an evolution from /i p to /j,p, with a NLO 

,2 /,,2 \ j f„ l„„j;v,„ /ATT T rjn- ln n - 1 f',,2 



kernel, allows one to resum leading (LL) a s hi n (fj, F / /j,q F ) and next-to-leading (NLL) agin 
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logarithms (collinear resummation) 
can be found, e.g., in Ref. [2]. 

In particular, in H — > bb processes, for an evolution from /x f — mb to fip ~ tuh, the large logarithms 
\n(m 2 H I 'ml) appearing in the massive calculation (2) are resummed with NLL accuracy. Likewise, since 
the perturbative fragmentation approach and the collinear resummation are process-independent, we can 
resum, following the same procedure, NLL logarithms \n(s/m\) and ln(raf/m^) in the massive b spectrum 
in e + e~ annihilation [2] and top decay [1]. 

Furthermore, both the coefficient functions in Eq. (5) and Refs. [1, 6] and the initial condition of the 
perturbative fragmentation function (6) present terms that become large once the 6-quark energy fraction 
Xb approaches 1, which corresponds to soft-gluon radiation. Soft contributions to the initial condition are 
process-independent and were resummed in [6] in the NLL approximation. NLL soft-gluon resummation 
in the coefficient function is instead process-dependent and it was implemented in e + e~ annihilation [6], 
top decay [9] and H — > bb [3] processes. For the sake of brevity, I do not report here the expression for 
the soft-resummed initial condition and coefficient functions. 

In Fig. 1, the bottom energy distribution in Higgs decay is given; the result is qualitatively similar to 
what can be found in e+e - processes or top decay. I have plotted the NLO massive calculation (dashed), 
including NLL collinear resummation (dotted), and with both NLL collinear and soft resummations 
(solid). I have set fi = fip = tuh, [i = A*of = m b, % = 120 GeV, mb — 5 GeV. The fixed-order 
calculation lies below the resummed ones and grows as Xb — > 1, since <~ 1/(1 — the collinear-resummed 
spectrum exhibits instead a sharp peak at large Xb- Soft resummation is relevant at Xb > 0.6: the 
distribution is further smoothed and shows the Sudakov peak at Xb — 0.97. References [1, 3, 6] show that, 
after collinear and soft logarithms are resummed, the Xb distribution exhibits very mild dependence on 
factorization and renormalization scales, which corresponds to a reduction of the theoretical uncertainty 
on the prediction. 

It is interesting to compare the b spectrum in the considered three processes. In Fig. 2 we plot the 
Xb distributions in H — > bb, e + e~ — > bb at ^/s = mz, i.e. LEP I and SLD ccntre-of-mass energy, and 
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Figure 1: b spectrum in Higgs decay according to the NLO massive calculation (dashed), including NLL 
collinear resummation (dotted), and both NLL collinear and soft resummations (solid). In the inset, the 
same curves are shown at large Xb, on a logarithmic scale. 




Figure 2: b spectrum in Higgs decay, e + e annihilation at = 91.2 GeV (dashes), top decay (dots). 
All curves include NLL collinear and soft resummation. 
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t — > bW. The b spectrum in H decay is the highest at small Xb and the lowest at large Xb\ in top decay 
it is shifted toward large x\, and sharply peaked very close to 1. The e + e~ — ► bb prediction is similar 
to the one yielded by Higgs decay, although some difference is still visible around the Sudakov peak, ft 
was shown in [3] that setting yfs = tuh makes the predictions in H — > bb and e + e~ annihilation more 
similar still. The difference noticed between the b spectra in Higgs and top decay can be experimentally 
relevant when distinguishing jets initiated by b quarks in the process pp — > ttH, followed by H — > bb. In 
fact, this is one of the search channels of the Standard Model Higgs boson at the LHC, for a Higgs mass 
m H < 135 GeV [10]. 

I would like to present results for 6-fiavoured B hadron production in the investigated processes. 
The B spectrum will be given by the convolution of the 6-energy distribution with a non-perturbative 
fragmentation function associated with the hadronization. As in [1, 3, 9], I consider a power law with 
two tunable parameters 

Dnp (^) = BW+ l a+1) (i-*r^ (7) 

the model of Kartvelishvili et al. [11] 

D n P(x;6) = {l + S)(2 + d)(l-x)x s , (8) 

and the Peterson model [12] 

D np & c) - n 71 A m — w ( 9 ) 

x[l — 1/x — e/(l — x)Y 

In Eq. (7), B(x, y) is the Euler beta function; in (9) A is a normalization constant. 

The parameters a, (i, 5 and e can be obtained after fitting models (7)-(9) to xb data in e+e~ 
collisions, where ig is the normalized B energy fraction. In [3, 15], data from ALEPH [13] and SLD [14] 
Collaborations were considered. In [15] it was reported that if we fit independently the non-perturbative 
models to ALEPH and SLD, we get pretty different best-fit parameters. Unlike ALEPH, which detected 
B mesons, SLD was able to reconstruct some B baryons as well; however it is a very small fraction of 
the whole sample and that may not be the reason for the discrepancies discussed in [15]. 

In [3] a combined fit to ALEPH and SLD was performed instead, as if all data came from the same 
sample. This approach is used here too. In the fits, the correlations among data points are neglected, 
and the considered data are in the range 0.18 <x B < 0.94. The best-fit parameters, along with the x 2 
per degree of freedom, are listed in Table 1. 



a 


0.90 ±0.15 


a 


16.23 ± 1.37 


X >,/3)/dof 


33.42/31 


<5 


17.07 ±0.39 


xWdof 


33.80/32 


e 


(1.71 ±0.09) x 10- a 


rWdof 


166.36/32 



Table 1: Results of combined fits to e + e~ — > bb data from the ALEPH and SLD Collaborations, using NLO 
coefficient functions, NLL DGLAP evolution and NLL soft-gluon resummation. I have set A = 200 MeV, 
Mof = Mo = m b = 5 GeV and fip = fi = = 91.2 GeV. 

The power law (7) and the Kartvelishvili model (8) yield very good fits, while the Peterson model (9) 
is unable to reproduce the data. In [3] it was also shown that the inclusion of NLL soft resummation in 
the coefficient function and in the initial condition of the perturbative fragmentation function is necessary 
to reproduce the data. Figure 3 shows the data samples, along with the fitting curves corresponding to 
the central values of a, (i, 5 and e quoted in Table 1. We see that models (7) and (8) yield approximately 
the same distribution: in fact, from Table 1 we learn that, within the error ranges, a is consistent with 
1 and P with 5. The Peterson non-perturbative fragmentation function is unable to fit the data. 
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Figure 3: ALEPH and SLD data on 5-hadrons, along with the best-fit curves according to the power 
law (7) (solid line), the Kartvelishvili model (dashes) and the Peterson model (dots), with a, (3, 5 and e 
given by the central values reported in Table 1. All curves include NLL collinear and soft resummation 
in the parton-level calculation. 

Given the results of Table 1, we can predict the i?-hadron spectrum in Higgs or top decay. In Fig. 4 
we show the xb distribution in Higgs decay according to models (7) and (8). Plotted are the edges of 
bands at one-standard-deviation confidence level. We can see that the predictions yielded by the two 
models are in statistical agreement. 

In Fig. 5 we compare 6-flavoured hadron spectra in e + e~ annihilation at LEP energy, Higgs and top 
decay, according to the power law (7). The results are similar to the ones at parton level presented in 
Fig. 2: the top-decay spectrum is shifted toward larger xb, the H — > bb one is the highest at low xb, the 
e + e~-prediction annihilation lies between the two. 

Finally, we can use the moments on B production in e + e~ annihilation from the DELPHI Collab- 
oration [16] to predict the same moments in top and Higgs decay. In iV-space we have: = o h n D n ^ , 
Tfj = r b N D^ = r^CT^/er^r, where <jn and are the moments of e + e~-annihilation cross section, Higgs 
or top width, at hadron or parton level. The results are shown in Table 2: the moments of the top-decay 
width are the largest, the ones of the Higgs width are the lowest, and the e + e~-annihilation ones are in 
the middle. This result is consistent with what is observed in x^-space in Fig. 5. 
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(ar*) 


<* 4 } 


e + e data 


0.7153±0.0052 


0.5401±0.0064 


0.4236±0.0065 


0.3406±0.0064 


e+e- NLL a b N 


0.7801 


0.6436 


0.5479 


0.4755 


D J [B] 


0.9169 


0.8392 


0.7731 


0.7163 


t-decay NLL T b N 


0.7884 


0.6617 


0.5737 


0.5072 


t-decay 


0.7228 


0.5553 


0.4435 


0.3633 


iJ-dccay NLL T b N 


0.7578 


0.6162 


0.5193 


0.4473 


i7-dccay T% 


0.6948 


0.5171 


0.4015 


0.3204 



Table 2: Experimental data for the moments from DELPHI [16], the resummed e + e~ perturbative calculations 
for a% [6], the extracted non-perturbative contribution D^f . Using the resummed perturbative result on top 
and Higgs partonic widths, a prediction for the moments Ff, t — » bW and H — > bb processes is given. 

In summary, I have reviewed recent results on bottom-quark fragmentation in e + e~ annihilation, 
Higgs and top decay. Parton-level spectra exhibit a relevant impact of the inclusion of NLL collinear 
and soft resummation; using information from LEP and SLD experiments, ^-flavoured hadron energy 
distributions have been predicted. 
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Figure 4: £>-hadron spectrum in Higgs decay: the hadronization is described according to the power law 
(7) (solid) and the Kartvelishvili model (8) (dashes). Plotted are one-standard-deviation bands for the 
non-perturbative parameters a, [3 and 5. 
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